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ABSTRACT 

Context. Prominences are partially ionized, magnetized plasmas embedded in the solar corona. 
Damped oscillations and propagating waves are commonly observed. These oscillations have 
been interpreted in terms of magnetohydrodynamic (MHD) waves. Ion-neutral collisions and 
non-adiabatic effects (radiation losses and thermal conduction) have been proposed as damping 
mechanisms. 

Aims. We study the effect of the presence of helium on the time damping of non-adiabatic MHD 
waves in a plasma composed by electrons, protons, neutral hydrogen, neutral helium (Hei), and 
singly ionized helium (He n) in the single-fluid approximation. 

Methods. The dispersion relation of linear non-adiabatic MHD waves in a homogeneous, un- 
bounded, and partially ionzed prominence medium is derived. The period and the damping time 
of Alfven, slow, fast, and thermal waves are computed. A parametric study of the ratio of the 
damping time to the period with respect to the helium abundance is performed. 
Results. The efficiency of ion-neutral collisions as well as thermal conduction is increased by 
the presence of helium. However, if realistic abundances of helium in prominences (~ 10%) are 
considered, this effect has a minor inffuence on the wave damping. 

Conclusions. The presence of helium can be safely neglected in studies of MHD waves in par- 
tially ionized prominence plasmas. 

Key words. Sun: oscillations - Sun: magnetic fields - Sun: corona - Sun: prominences 

1. Introduction 

Small-amplitude oscillations and propagating waves are commonly observed in both quiescent 
and active region prominences/filaments. They have been interpreted in terms of magnetohydrody- 
namic (MHD) eigenmodes of the magnetic structure and/or propagating MHD waves. The reader is 
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referred to some recent reviews for more information about the observational and theoretical back- 



grounds dOliver & Ballestei 



2002; 



2008) 



En g void 



2004 



Ballestei 



2006; 



Baneriee et al 



2007 



En g void 



Prominence oscillations are known to be quickly damped, with damping times correspondin g 



Ohver 



2009; 



Mackav et al 



2009h . 



to a few oscillatory periods (this topic has been reviewed by 
Several damping mechanisms of MHD waves have been proposed, non-adiabatic effects and ion- 
neutral collisions being the more extensively investigated. In order to understand in detail these ef- 
fects, they have been s tudied in simple configurations such as unbounded and homogeneous media. 



Carbonell et al 



( 120041) investigated the time damping in a homogeneous prominence medium tak- 
ing non-adiabatic effects (optically thin radiation losses and the rmal conduction) into account. Later 
on, the spatial dampi ng was studied by Carbonell et alj (120061) and the effect of a background mass 



Carbonell et al 



(120091) . Subsequently, some works have extended these pre - 



flow was analyzed by 
vious results by considering the presence of the coronal medium (ISoler et al. 



2007 



2008 



2009ah . 



The common conclusion of these investigations is that only slow and thermal waves are efficiently 
damped by non-adiabatic effects, while fast waves are very slightly damped and Alfven waves are 
completely unaffected. 

On the other hand, the influence of partial ionization on the propagation and tim e damping of 
MHD waves has been also inve stigated in an unbounded medium. 



Forteza et al 



(|2007|) followed the 



treatment bv lBraginskiil d 1965b and derived the full set of MHD equations al ong with the dispe rsion 



Pinto et al 



2008). The 



relation of linear waves in a partially ionized, single-fluid plasma (see also[ 
presence of electrons, protons, and neutral hydrogen atoms was taken into account, whereas helium 
and other species were not considered. In a subsequent work dForteza et alj|2008l) . they extended 



their previous analysis by considering radiative losses and thermal conduction by electro ns and 



neutrals. Their main results with respect to the fully ionized case dCarbonell et al 



2004) were. 



first, that ion-neutral collisions (by means of the so-called Cowling's diffusion) can damp both 
Alfven and fast waves but non-adiabatic effects remain only important for the damping of slow 
and thermal waves, and second, that there exist critical values of the wavenumber in which the 
real part of the frequency vanishes, so wave propagation is not possible for larger wavenumbers 



Again, a pplications to a more complex cylindrical geometry have been also performed dSoler et al 
2009blcl) 



On the basis of these previous results, it seems clear that partial ionization plays a relevant role 
on wave propagation in prominences. Prominences are roughly composed by 90% hydrogen and 
10% heUum but, to date, all the investigations considered a pure hydrogen plasma. Therefore, the 
effect of the presence of helium on the propagation and damping of MHD waves is still unknown 
and is the motivation for the present work. Here, we consider an unbounded and homogeneous 
prominence medium permeated by a homogeneous magnetic field. The plasma is assumed to be 
partially ionized, electrons, protons, neutral hydrogen, neutral helium ( Hei), and singly ionized 



Gouttebroze & Labrosse 



helium (Hen) being the species taken into account. Recent studies by 
d2009l) indicate that for central prominence temperatures, the ratio of the number densities of He n 
to He I is around 10%, whereas the presenc e of He m is n e gligible. Thi s result allows us to neglect 
He m in this work. Extending the works by 



Forteza et al 



d2007 



2008), the derivation of the basic 



MHD equations for a non-adiabatic, partially ionized, single-fluid plasma has been generalized by 
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considering now five different species, allowing us to study how the presence of neutral and singly 
ionized helium affects their previous results. 

This paper is organized as follows. The description of the equilibrium and the basic equations 
are given in Sect.|2] The results are discussed in Sect. [3] Finally, Sect. |4]contains the conclusion of 
this work. 



2. Equilibrium and basic equations 



Our equilibrium configuration is a homogeneous and unbounded partially ionized plasma com- 
posed by electrons, protons, neutral hydrogen, neutral helium, and singly ionized helium. Hereafter, 
subscripts e, p, H, He i, and He ii explicitly denote these species, respectively. The magnetic field 
is also homogeneous and orientated along the jc-direction, Bq — BqC x, with Bq - 5 G. We adopt 
the single-fluid approximation. Following 



Forteza et al 



(120071 



20081) and neglecting the electron 



contribution, we define the center of mass velocity, v, as follows. 



V ~ ^pVp + ^hVh + ^HelVHel + ^HellVHell, 



(1) 



with the relative density of species a, and v„ the corresponding species velocity. Equivalently, 
the equihbrium total density, po, and gas pressure, pQ, are. 



PO ~ Pp +Ph +PHei +PHen, 



(2) 
(3) 



Since pa - ^aPa, we get the relation + + ^Hei + ^Hen ~ 1- We assume a strong thermal 
coupling between species, so all the species have the same equilibrium temperature Tq. Then, the 
three equilibrium quantities are related as follows. 



R 

Po = Po — Tq, 

where R is the ideal gas constant and fi is the mean atomic weight, 
1 



2fp +^H + J^Hel + 5^HeII 

With the help of some definitions, we can express pin a more convenient form, 
Ah 



H = 



1 - [(1 + ^He) + (1 + 25He) jPh] ^Hel ' 



with 



J"H = 



^p+^H 



^H 



(4) 



(5) 



(6) 



(7) 



2^p + ^H ^He I 

The quantity fiu is equiva l ent to the mean atomic weight of a pure hydrogen plasma defined 



in Eq. (3) of 
plasma and fin 



Forteza et al 



(120071) . and ranges between pH - 0.5 for a fully ionized hydrogen 
1 for a fully neutral hydrogen gas. On the other hand, ^He indicates the he- 
li um ionization degree. A reahstic value of this parameter is ^He = 0.1 according to the results 
of iGouttebroze & Labrossd ( 120091) . From Eq. (|6j one can see that p > pn due to the presence of 



helium. In the absence of helium, ^Hei - ^Heii = so p - pu. Figure [T] displays the dependence of 
p on ^Hei for several values of pu. 



4 
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Fig. 1. Mean atomic weight, jl, as a function of the relative neutral helium density, ^He i, for flu - 0.5 
(dotted line), yUe = 0.6 (dashed line), /ie = 0.8 (solid line), and /ie = 0.95 (dash-dotted line). In all 
cases, ^He = 0.1. 



The details of the derivation of the basic gov erning equations f or a non-adiabatic, partiall y 



ionized, one-fluid plasma can be followed in, e.g.. lBraginskiil (Il965b : 



Pinto et al. 



Forteza et al 



l atic, p : 
( 20071 



(l2008b . Here, we follow the same procedure but generalize the analysis of lForteza et al 



2008) 



( l2007b by including additional species. In brief, the separate governing equations for the five 



species are adde d and a generaUzed Ohm's law is obtained. These basic equations correspond 
to Eqs. (l)-(6) of lForteza et al.l (l2008h . which are formally identical in our case. A key step in the 



present derivation is to compute the density current, j, as 



(8) 



along with the condition = "p + "Hen, where n^, rip, and neeii are the electron, proton, and 
He 11 number de nsities, respectively, and e is the electron charge. The resulting induction equation 
(see Eq. [14] of 



Forteza et al. 



2007h contains several diffusion terms whose coefficients depend 



on the collisional fre quencies between species. The physical meaning of these nonideal terms is 
explained in detail in 



Pinto et al 



(I2008h . In particular, ion-neutral collisions are responsible for the 
so-called Cowling's diffusion, which is much more efficient than Ohm's diffusion in a partially 
ionized plasma. However, some terms are not relevant f or our present application. Hall's effect is 
negligible in prominence conditions dSoler et al.ll2009d) . and the so-called "Biermann's battery" 
term is identically zer o in a homogeneous medium. For this reason, our final form of the induction 
equation (Eq. [21] of lForteza et al.luOOTb only contains the terms corresponding to Ohm's and 



ambipolar (Cowling's) diffusion, along with the diamagnetic current term. 

Ohm's, T], and Cowling's, rjc, coefficients of magnetic diffusion can be expressed in terms of 
their corresponding conductivities. 



1 

'7 = — > 



1 



(9) 
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with /i = 471 X 10"^ N A"^. So, Ohm's and Cowling's conductivities, as well as the diamagnetic 
current coefficient, H, which applies in our case when helium is included are: 

cr 



2 2 

e nt 



„ (^H+^Hel) 



0"C 



//an 



^p^H - ^^H^Hell + ^^p^Hel + ^^Hel^Hell 



(10) 



(11) 



In addition, a^, a^n, and are the electron, electron-neutral, and neutral friction coefficients, re- 
spectively, whose expressions depend on the sum of the friction coefficients between particular 
species. 



ae = ^ep + O'eH + CeHel + CeHell, 

an - O'eH + aeHel + ^pH + ttpHe I + ffHe IlH + aHellHel- 

Each particular friction coefficient, a^/j', is computed as 



(12) 
(13) 
(14) 

(15) 



with tip the number density of the species fi, vpp' the coUisional frequency between species /3 and 
/?', and 

nipnip' 



De Pontieu et al 



nip + nipi 

with ffifl the mass particle of the species /3. As given by 
(l2009bl) . the coUisional frequencies between electrons and protons or He ii are 

Vei = 3.7 X 10" 



(16) 



( l200lb . see also 



Soler et al 



_gni In A 

7^3/2 ' 



(17) 



with i = p or He II and In A the Coulomb logarithm, while the coUisional frequency between a 
charged species, q = e, p, or He ii, and a neutral species, n = H or He i, is 



Vqn = «n 



(18) 



''qn 



with the Boltzmann's constant, and Sqn the coUisional cross-section. Here, we consider the 
values Sen - 10"'^ m^, and Eon = ^He 



5 X 10-''' m^. 



On the other hand, the thermal conductivity due t o neutra l s (Eq. [16] of lForteza et alj|2008l) has 
to include now the helium contribution. According to lParken (Il953h . a corrected expression for the 
conductivity of neutrals in MKS units is 



= + ^Hel = (2.44 X lO^^^H + 3.18 X lO-^^Hel) T^^l 



(19) 



Finally, we assume a n optically thin radiati on dHildneiil 19741) to represent the hydrogen radiative 



losses. According to ICox & Tucked (119691 see their Fig. 3), the radiative losses by helium are 
several orders of magnitude smaller than those of hydrogen for typical prominence temperatures 
(~ lO'* K) and, therefore, irrelevant for the present invest i gation . 



Hereafter, our analysis follows that of 



Forteza et al. 



( 20081) . We linearize the basic equations 



and assume small perturbations proportional to exp (icot + ik^x + ik,z). Then the resulting equations 



6 
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Fig. 2. Ratio of the damping time to the period, tu/P, versus the wavenumber, k, corresponding 
to the (a) Alfven wave, (b) fast wave, and (c) slow wave for - tt/A, fin = 0.8, and ^He =0.1. 
The different linestyles represent ^Hei - 0% (solid line), ^Hei = 10% (dotted line), and ^Hei = 20% 
(dashed line). The results for ^Hei = 10% and ^He = 0.5 are plotted by means of symbols for 
comparison. The shaded regions correspond to the range of typically observed wavelengths of 
prominence oscillations. 



(Eq. [18]-[27] of lForteza et al.ll2008h are combined and finally two difi^erent, uncoupled dispersion 
relations, one for Alfven waves (their Eq. [28]) and another for magnetoacoustic and thermal waves 
(their Eq. [30]) are obtained. Note that although our definitions of 77, 77c, H, and Kg contain the 
effect of helium, the resulting dispersion relations are formally identical to those of lForteza et al 
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( l2008h . For the sake of simp licity, we do not write again these expressions here and refer the 



Forteza et al 



(|2008|) . The dispersion relations are numerically solved for real values of 



reader to 

the wavenumber modulus, k - + kl, and the angle 6 between Bq and k. A complex frequency, 
CD - (Wr+Zwi, is obtained. The period, P, and damping time, td, are related to the real and imaginary 
parts of the frequency as follows, 

2n 1 

P=—, TD = — . (20) 

3. Results 

In the following computations, we consider typical prominence conditions, po = 5 x 10"" kg m^^* 
and Tq - 8000 K. Quantities /jh, ^hci, and ^He are considered free parameters. We focus our 
attention on the effect of the relative neutral helium density, ^Hei, on the ratio td/P. 

3.1. Free propagation in an unbounded medium 

First, we assume 9 - jijA. Figure|2]displays t^IP as a function of k for the Alfven, fast, and slow 
waves. The results corresponding to several helium abundances are compared for hydrogen and 
helium ionization degrees of /in — 0.8 and ^He =0.1, respectively. We see that even in the case of 
the largest quantity of helium considered (^He i = 20%), the presence of helium has a minor effect 
on the results. In the case of Alfven and fast waves (Fig.|2^,b), their critical wavenumber (i.e., the 
value of k which causes the real part of the frequency to vanish) is shifted toward slightly smaller 
values. So, the larger ^Hei, the smaller k^. This result ca n be understood by c onsidering that the 
Alfven wave critical wavenumber, k^, given by Eq. (38) of 



Forteza et al. 



mm is 



= (21) 
{tjc + rj tan^ 9) cos 9 ' 

with VA = Bo/ yfjjpo the Alfven speed. Equation (1211 1 is also approximately valid for the fast wave 
critical wavenumber Then, we see that k^ is inversely proportional to Cowling's diffusivity, rjc. 
Since rjc is larger in the presence of helium than in the pure hydrogen case due to additional 
collisions of neutral and singly ionized helium species, k^ is therefore smaller. Turning our atten- 
tion to the slow wave (Fig. |2};), we see that the maximum and the ri ght-hand side minimu m of 
Tp/P are also slight ly shifted toward smaller values of k. Results from 



Carbonell et al 



Forteza et al. 



(|2004 and 



(120081) indicate that thermal conduction is responsible for these maximum and min- 
imum of td/P. Thus, the additional contribution of neutral helium atoms to thermal conduction 
(Eq.[T9Tl causes this displacement of the curve of tu/P. As for Alfven and fast waves, this effect is 
of minor importance. For comparison, equivalent results with ^He i = 10% and ^He = 0.5 are plotted 
by means of symbols in Fig.|2] We see that for realistic values of ^He, its role is almost irrelevant, 
meaning that the presence of Hen can be neglected. It is worth mentioning that we have repeated 
these calculations for other values of fiu and similar results have been obtained. 

Next, we study the thermal mode. Since it is a purely damped, non-propagating disturbance 
(wr - 0), we only plot the damping time, td, as a function of k for jlu = 0.8 and due =0.1 
(Fig. [3j. We can see that the effect of helium is different in two ranges of k. For k > 10 m"', 
thermal conduction is the dominant damping mechanism. So, the larger the amount of helium, 
the smaller tq because of the enhanced thermal conduction by neutral helium atoms. On the other 
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hand, radiative losses are more relevant for k < 10 m ' . In this region, the thermal mode damping 
time grows as the helium abundance increases. Since these variations of the damping time are very 
small, we have to conclude again that the damping time obtained in the absence of helium does not 
significantly change when helium is taken into account. Computations with other values of jlu and 
^He do not modify this statement. 
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Fig. 3. Damping time, td, of the thermal wave versus the wavenumber, k, with = 7r/4. The 
diff'erent hnestyles represent: ^Hei = 0% (solid line), ^Hei = 10% (dotted line), and ^Hei - 20% 
(dashed line). In all computations, /zh = 0.8 and ^hc = 0.1. The result for ^Hei - 10%and5He = 0.5 
is plotted by means of symbols for comparison. 



3.2. Constrained propagation by a waveguide 



We can estimate the effect of a magnetic structure, say a slab or a cylinder, which would act as a 
waveguide. To do so, we set the wavenumber component in the perpendicular direction to mag- 
netic field fines to a fixed value, k^L - jt/2, with L a typical lengt h-scale in the perpend icular 



Lin et al 



2007 



20091) show 



direction. Since high-resolution observations of filaments (see, e.g., 
fine-structures (threads) with a typical width of ~ 100 km, we select L = 10^ m as our perpendicu 
lar length-scale. Therefore, the propagation angle 6 depends now on kx. 



6 = arctan 



k,L 



(22) 



Figure H] displays the results for the Alfven, fast, and slow waves. We see that the behavior of 
the three solutions is substantially diff'erent from that of the free propagation case. The Alfven 
mode now possesses an additional critical wavenumber for small values of kx, namely k'^T , which 
is indepe ndent of the ioniza tion degree and the helium abundance. It can be approximated as (see 
Soler et al.ll2009b[ Eq. [38]), 



details in 



fe? = 

2va ~ 8vAi 



(23) 



On the other hand, the fast wave is now more attenuated in the relevant range of wavenumbers than 
in the free propagation case, whereas the slow wave also has a new critical wavenumber, namely 
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Fig. 4. Ratio of the damping time to the period, tu/P, versus the wavenumber component parallel 
to magnetic field lines, fc^, corresponding to the (a) Alfven wave, (b) fast wave, and (c) slow wave 
for k,L = 7r/2, with L = 10^ m, jlu - 0.8, and ^He = 0.1. The diff'erent linestyles represent ^Hei = 
(solid line), ^Hei = 10% (dotted line), and ^Hei - 20% (dashed fine). The vertical dot-dashed lines 
in (fl) and (c) correspond to the approximated critical wavenumbers given by Eqs. ( l23T l and ( |241) . 
respectively, for ^He I = 10%. 



kf which falls within the re levant range . An exp ression for the slow mode critical wavenumber is 



also provided by Eq. (48) of 

^2 



Soler et al 



(l2009bl ). which in our present notation is. 



2vl ' 



8viL2 



(24) 
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where Cs = -^Jypo/po is the sound speed, with y = 5/3 the adiabatic index. The slow mode critical 
wavenumber is shifted toward larger values as the heUum abundance increases. Note that there is 
no additional critical wavenumber for the fast wave. These approximated critical wavenumbers are 
indicated by means of vertical lines in Fig. |4] We see an excellent agreement in the case of the 
Alfven mode critical wavenumber (Eq. ||23l ). For the slow wave, the approximated value (Eq. Il24l ) 
is slightly larger than that numerically obtained. 

Finally, we have also computed the results in the case of the guided thermal disturbance. We 
find that the thermal mode behavior is the same in the waveguide case and in the free propagation 
case. Hence, this mode is not affected by the variation of the propagation angle and no further 
comments are needed. 



4. Conclusion 

In this work, we have studied the effect of helium (He i and He ii) on the time damping of thermal 
and MHD waves in a part i ally ionized prominence plasma. This is an extension of previous inves- 



tigations bv lForteza et al.l ( 12007 



2008h in which helium was not taken into account. We conclude 



that, although the presence of neutral helium increases the efficiency of both ion-neutral collisions 
and thermal conduction, its effect is not important for realistic helium abundances in prominences. 
In addition, due to the very small He n abundance for central prominence temperatures, its presence 
is irrelevant to the wave behavior This conclusion applies both to the free propagation case and 
the constrained propagation by a waveguide case. Although the role of He ii (or even He iii) could 
be larger for typical prominence-corona transition region temperatures, the present result allows 
future studies of MHD waves and oscillations in prominences to neglect the presence of helium. 
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